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Abstract 

We consider multidimensional attractive particle systems with prod- 
uct invariant measures and open boundaries. We show that the Euler 
hydrodynamic limit is a scalar conservation law with boundary condi- 
tions of Bardos, Leroux and Nedelec. The stationary profile between 
parallel hyperplanes follows by proving asymptotic stability of entropy 
solutions. Examples include fluxes with several minima and maxima. 
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1 Introduction 

Stochastic lattice gases in contact with reservoirs are relatively tractable 
and thus widely studied examples of nonequilibrium stationary states. The 
derivation of the stationary macroscopic profile (hydrostatic limit) is a natu- 
ral problem in this context. For diffusive bulk dynamics, robust methods (see 
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e.g. [221 12SI E3 SHI) show that the hydrostatic profile is the stationary solu- 
tion to the hydro dynamic equation (a possibly nonlinear diffusion equation) 
with Dirichlet boundary conditions imposed by the reservoir densities. For 
driven lattice gases, no general approach is available. The one-dimensional 
asymmetric exclusion process was addressed in [16] and [121 H6] by model- 
specific methods. The hydrostatic profile consists of three uniform phases 
outside the phase transition line A; < A n A/ + A r = 1, where A/ and A r denote 
the density of the left and right reservoir. On the phase transition line it 
consists of a randomly located shock connecting A; and A r . In both cases 
this profile is a stationary solution to the hydrodynamic equation, now of 
hyperbolic type: 

d t p(t,x) + V x .f(p(t,x))=0 (1) 

with the current- density function given by f(p) = p(l — p). However, out- 
side phase transition, it cannot satisfy Dirichlet boundary conditions at both 
ends. In this paper we provide the proper interpretation of the hydrostatic 
profile in terms of boundary conditions and use this new point of view to rig- 
orously extend the results of [JB] and [HI HS| to more general driven particle 
systems and higher dimensions. In [53], a variational formula was introduced 
heuristically for the bulk density of one-dimensional driven lattice systems 
with open boundaries, governed by a hydrodynamic equation of the form 
([1]). When the current f(p) exhibits two local maxima with a local minimum 
inbetween (as in the KLS model, see [38J), this gives rise to a seven-phase 
diagram. There is so far no mathematical proof of this result, nor any at- 
tempt to generalize it to higher space dimension. We establish such results 
for a wide class of models as a consequence of the hydrodynamic limit. 

The hydrodynamic behavior of particle systems with open boudaries is so 
far understood in diffusive regimes ([231 SHI E]), but not in Euler scaling, due 
to the specific nature of the boundary conditions for ([1]). We prove here that 
the Euler hydrodynamic limit is the entropy solution to an initial-boundary 
problem for (fl|) with the BLN boundary conditions introduced in [7], that 
are relevant for hyperbolic conservation laws. Instead of fixing the bound- 
ary value like Dirichlet conditions, these conditions impose a set of possible 
boundary values depending on the boundary datum ([IH]). Our result is 
established for attractive particle systems on Z d with product invariant mea- 
sures, that is in framework of the celebrated paper [55], which proved that 
the hydrodynamic limit on the whole space was the entropy solution to ([!]). 
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Here we restrict particles to the lattice discretization of a regular open sub- 
set Q of M. d . Boundary dynamics are defined to model exchanges of particles 
with a reservoir of given density profile A(.) G L°°(dQ). Our proof uses 
a generalized formulation Q64J) of the BLN boundary conditions that does 
not explicitely involve a trace for the solution. Producing this formulation 
from the microscopic boundary dynamics involves adequate coupling of open 
systems. We next derive the hydrostatic profile and local equilibrium in a do- 
main Q between two parallel hyperplanes coupled to uniform reservoirs. The 
result is a <i-dimensional version of [53] , and holds more generally between 
perturbed hyperplanes away from the perturbation. Other geometries could 
be tractable by our method and induce a new type of phase diagrams: we 
hope to obtain rigorous results in further work ([S])- The hydrostatic limit is 
derived from the hydrodynamic limit after proving a large-time convergence 
theorem for entropy solutions with BLN boundary conditions. The latter 
problem has been studied only for fluxes that are convex ([49]) or with a 
single maximum ([50]). By comparison and entropy dissipation arguments 
we treat general fluxes. 

It would be interesting to treat other models with a hydrodynamic limit 
of the form (pQ). These include systems with a microscopic Lax-Hopf for- 
mula, e.g. the totally asymmetric i^-exclusion pro cess ([58]); attractive sys- 
tems with non-explicit invariant measures ([SS],[S]); KLS and other nearest- 
particle models ( [M] , [E] , [2] , [3] ) ; weakly asymmetric systems with vanish- 
ing viscosity (|25j-[2H]). Boundary extension of the Lax-Hopf formula ([35]) 
could be sought for i^-exclusion process and related models. For the models 
in [58| 156} [5], invariant measures are not explicit, and not even known to 
exist for all densities. One is led to define simple boundary mechanisms that 
are not related to invariant measures. In |3T] this was conjectured to induce 
a non-explicit effective boundary density. We believe that the hydrodynamic 
limit would still given by BLN conditions for this effective boundary density. 

The paper is organized as follows. In section [2] we define the framework, 
state results, and give examples. The hydrodynamic limit (Theorem 12. ip is 
established in Section [31 the hydrostatic limit (Theorem 12. 3j) in Section [U 
and large-time convergence of entropy solutions (Theorem 12.21) in Section [5j 
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2 The framework and results 



Notation. Let N = {1, 2, . . .} and Z+ = {0, 1, . . .}. For x, y G R d and 8 > 0, 
|x| = maxj = i r .. ) d \xi\, x A y = mm(x,y), x V y — max(x, y), x + = x V 0, 
x~ = -(x A 0), B 5 (x) = {z G R d : |z - x| < 5}. 7i d_1 denotes the (d - 1)- 
dimensional Hausdorff measure in M. d . E := E z is the set of particle configu- 
rations rj = (t](x) : x G Z d ), where = [0, /C] flZ + , and /C G NU {+00} is the 
maximum number of particles per site. A local (or cylinder) function /onE 
is a function that depends only on the restriction of 77 to a finite subset of Z d . 
The partial product order on E, i.e. rj < £ iff. rj(x) < for all x G Z d , in- 
duces (see e.g. [SHI SI] ) a partial stochastic order among probability measures 
on E: for two probability measures \L\ and on E, \i\ < /x 2 iff-, equivalently: 
(a) For every nondecreasing function / on X, f fdfii < f fdfi 2 ', (b) There 
exists a probability measure fj, on X 2 , with marginals fix and 112, supported 
on {(771,772) G E 2 : T]\< r] 2 }. Spatial shift on E, defined by r x r]{.) = r](x + .) 
for x G Z d , is extended to functions / G M E by r x f = f o t x , and to operators 
L on R E by t x L = L o r x . 

Framework. We consider particle systems on Z d (i.e. Feller processes on 
E, with generator denoted by L and semigroup denoted by (e tL ,t > 0)) 
with the following properties - see examples in Subsection 12. 21 (a) local and 
translation invariant dynamics, i.e. 

L = L x , L x = t x Lqt^ x (2) 

x&Z d 

where Lq is a bounded Markov generator acting only on B r (0) n Z d for some 
r G N (boundedness is not essential and could be replaced by a linear growth 
assumption on jump rates as in [1]); b) conservation of particles: 

Lo Yl V(x) = (3) 

a;GZ d : \x\<r 

with a bounded number of displaced particles, see (E2D ; (c) existence of a 
local monotone coupling (which implies attractiveness, see [B]); (d) an ex- 
ponential family of product invariant measures (z/ p ) pg 7^ on E, indexed by the 
mean particle density, p = ^[77(0)] G 1Z := [0, /C] PI [0, +00); by exponential 
we mean (this is in fact a consequence of ([3]) given (T5])), for 77, G Z + , 

1/^(0) =n) = u p M0)=n)ex P {[P(p)~P(p )]n-[A(P(p))-A(P(p ))]} 
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for p, p G 7Z\{0, fC}, where /?(.) is a C°°-diffeomorphism from 7Z\{0,JC} to 
(—00, /3 ); (e) irreducibility assumptions (see Subsection 12.21) . Let now par- 
ticles evolve in S C Z d . On the reservoir Z* d \S, we have a microscopic 
profile A(.), that is a bounded 1Z- valued field. We construct a Markov gener- 
ator Lg } x(.) 011 E s as follows. We decompose a particle configuration on Z d , 
77 G E, as 77 = f]° © where 77° G -E 5 is its restriction to S, and 77 G E z \ s its 
restriction to 7j d \S. If 77 is a E-valued r.v. with distribution p, p° denotes the 
distribution of 1] . p°, u°, etc., also denote probability distributions defined 
directly on E s . Let the reservoir state be the product measure t^s,\(.) on 
E zd \ s such that 

Vsx.)({rj(x) = n}) = 6 x{x) (n), \Jx G Z d \S, \Jn G Z+ (4) 
We set, for every local function r]° \— > f(rj°) on E s , 



Ls,x(.)f{r)°) = J WW® V)<&s,x(.)(v) 


(5) 


where Lf depends on 77 because of x near the boundary in 


©. © implies 


= E L - + E ^ 


(6) 


x£S r ° x€dS r 




where generators Lx ^ contain bulk-reservoir interactions. 


and, for 5 G N, 


S So := {xeZ d : d{x,Z d \S) > 5} 


(7) 


dS 5 := {xeZ d : d(x, S) < 5, d(x, Z d \S) 


< 5} (8) 


dS s+ := {xeZ d : d{x,S) < 5} 


(9) 


d(x, S) := inf{ \x — y : y G S} 


(10) 



A natural property of (J5j) is that a system in equilibrium remains so when in 
contact with a reservoir of equal density: 

Proposition 2.1 If A(.) = p, v° is invariant for Lg,x{.)- 

Indeed: by (JSJ) and invariance of u„ for L, for a local function /(77 ) on E s , 



L S)K) f(r,°)dv°M 



Lf(r)° ®rf)dv p {rf ©77) =0 
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2.1 Results 

Hydrodynamic limit. Let S = Qn '■= {x E 7* d : x/N £ fl}, where fl is 
an open subset of M. d with a boundary of locally finite perimeter ([H]), and 
N £ N a scaling parameter. Set E^ := E Qn and E N = E zd \ QN . We say 
a uniformly bounded sequence (Ajv(.), iV £ N) of microscopic reservoir pro- 
files on Z d \Q^ has limiting trace A(.) £ L°°(<9f2) (the macroscopic reservoir 
profile) on dfl, and we write Ajv(-) ~ A(.), if for some C > 0, 

limsupiV 1 ^ V ip(x/N)f(X N (x)) < C [ <p(x)f(\(x)) dH d - x {x) (11) 



for every nonnegative continuous functions / on 71 and <p with compact 
support on M. d . ffTTl) holds for instance if Aat(:e) = X(x/N), where A(.) a 
continuous function on R d \f2 with trace A(.) on dfl. To every configuration 
77 , we associate a sequence of empirical measures on fl: 

a N {r]°){dx) = N- 1 V°WjL(dx), N £ N (12) 

Let p No be, for each N £ N, a probability measure on E^. We say the 
sequence (p No )n has density profile p(.) £ L°°(0), if under p No the random 
measures a N (rj°)(dx) converge in probability as N — > oo to p(.)dx with re- 
spect to vague topology (i.e. against continuous test functions with compact 
support in fl). We write p N ° ~ p(.) if 

» N °(dri)= (g) vp {x) {drf(x)) (13) 

for some uniformly bounded family (p N (x), N £ N, x £ fijv) such that 

lim / \p N ([Nx\) - p{x)\dx = (14) 

N^oo J j 

for every bounded Borel set I G fl. Define 

(the microscopic and macroscopic flux), /i £ C 1 (7?-) because (z/ p ) is exponen- 
tial. Consider the scalar conservation law. 

d t p{t,x) + div x h(p(t,x)) = (16) 
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Theorem 2.1 Assume p No ~ p (.) G and \ N (.) ~ A(.) G L°°(<9fi). 

T/ien, /or every t > 0, p^t := p No e NtLn N- x N(-) has density profile pit, .), where 
p(., .) zs t/te unique entropy solution to [To]) in Q with initial datum po(.) an d 
boundary datum A(.). 

Following [61] , we say p(.) G L°°((0, +oo) x fi) is an entropy solution to (|T6l) 
in fi, with data p G and A G L°°(dQ), iff. there exists M > such 

that 



[<9 t i£>(i, x)(f)(p(t, x)) + V x (p{t, x).ip(p{t } x))] dtdx 

'(0,+oo)xR d 

+ Ml if{t,x)<f){X{x))dtdn d ~\x) + / (f{0,x)(p{p (x))dx > 

J(0,+oo)xdQ Jn 

(17) 

for every nonnegative smooth function with compact support in [0, +oo) x 
R d , and the Kruzkov entropy-flux pairs ([39]) (<f>,ip) = ((j>t,i>t) gi ven by 

<Pt{p) = (p- c) + , ij)+(p) = l ( o,+oo)(p - c)[h(p) - h(c)} - lg > 
<^(p) = (P - c )~; #T(p) = - 1 (-oo,o)(p - c)[h(p) - fc(c)] 



Existence and uniqueness of an entropy solution are proved in [64] . The 
original and equivalent definition of the entropy solution in [7] (see also [19] ) 
separately imposes (i) the classical Kruzkov's entropy conditions in Q: i.e., 
(ITTj) without the boundary term, for smooth ip with compact support in 
[0, +oo) x Q (here entropies C := 0+ + 0^7 would be sufficient); (ii) the BLN 
boundary condition 

p(t, x) G S\ { :\ n{x \ dt <g> H d -\dx)-&.e. on (0, +oo) x OVt (19) 

where n(x) is the inner unitary normal to <9f2, p(t, x) the trace of .), and 

fj[ (0 :={pGft: sgn(p-A)[/(p)-/(c)] < 0, Vc G [min(A, p), max(A, p)]} 

(20) 

is the set of admissible boundary values for a scalar flux function /(.) on 1Z 
and a density value A G TZ. Entropy conditions in Q ensure existence of a 
trace in sufficient sense to make ( Tl9l) meaningful (see [621 l4"0l 153]). 

Hydrostatic limit and stationary entropy solutions. Let n G R d be 

a unitary vector and a, 6 G R such that a < b. We assume now that f2 is a 
"perturbation" of 

fi£ 6 := {a; G R d : n.x G (a, 6)} (21) 



in the following sense: there exist — oo < a! < a <b <b' < +00 such that 

cnc qj )6 , (22) 

The boundary of Q is a disjoint union of components <9f2 7 for 7 G {a, 6}, such 
that n.x < a (resp. n.x > b) on dQ a (resp. on dQb)- If Q — fi™ b , <9f2 7 is the 
hyperplane {n.x = 7}. We consider ffl6|) in f2 with BLN boundary datum 

A(.) = A a lao a + (23) 

on <9f2, where (A a , A&) G 7£ 2 . For a scalar flux function / G C 1 (7?.), define 

A{ = inf{A < A a : / is constant on [A,A a ]} (24) 

X[ = sup{A > Xb '■ / is constant on [A&,A]} (25) 

if A a < A 6 , or 

A{ = sup{A > A a : / is constant on [A a ,A]} (26) 

Aj£ = inf{A < Xb : / is constant on [A, AJ} (27) 

if A a > Xb- For A a < Xb (resp. A a > Xb), let Ji4f(X a , Xb) denote the set of 
minimizers (resp. maximizers) of / on [A{ A X[, A{ V X[], and Rf(X a , Xb) its 
unique element when JHf(X a , Xb) is a singleton. 

Theorem 2.2 Assume h G C 1 ^), and Mh(.).n(X a , Xb) = {Rh(.).n(X a , Xb)}. 
Then, any entropy solution to [To]) in Q with boundary datum |13j) on dfl 
converges as t — ► 00 to Rh(\ n (X a , A ) in Ll oc (Q^ b ), and any stationary entropy 
solution p (.) G L°°(Q) is equal to Rh(.). n (X a , A ) in f2™ . In particular for 
Q = any entropy solution converges to the unique stationary solution 

Rh(.).n(X a , Xb). 

Remark, (i) The singleton condition is necessary for uniqueness of a sta- 
tionary solution in Q™ b . Otherwise, one can construct stationary solutions by 
connecting any sequence of increasing (decreasing) minimizers (maximizers) 
across hyperplanes parallel to the boundaries; in dimension one, these are 
the only stationary solutions, (ii) For h of the form (|T5|) in the framework 
of Section [21 h(.).n is not constant on any nontrivial interval unless it van- 
ishes on K. Then, given n G R d : either M h (.). n (X a , X b ) = {R h Q, n (X a , X b )}, 
A{ = A a and X[ = Xb for all (A a , A ); or Aih(.).n(X a , Xb) is a singleton for no 
(A a , Xb) (e.g. for the exclusion process with hyperplanes parallel to the drift, 
boundaries cannot determine inside density). 
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Theorem 2.3 Assume \ N {.) ~ A(.) G L°°(dQ), M h{ .), n (X a , X b ) = {R h Q. n (K, A 6 )}, 
and i^ is an invariant measures for Lfi Nt x N (.) satisfying 

There exists R G K and N G N suc/t toat > JV , v No < v° R (28) 

Then, as N — ► oo; (%) i/ie restriction to Q™b °f a: (77 ) defined by ( fJS|) 
/ias uniform density profile Rh(.).n{X a , Xb) (hydrostatic limit); (ii) for every 
bounded local function g on E, ~g N {x) := t^j^^^ )] converges in L 1 1 oc (f2^ 6 ) 
to <? := ^R fc( . ) .„(A a ,A 6 ) (»7°)] (foca* equilibrium). 

Remarks. (128]) is a microscopic counterpart of p s (.) G L°°(f2) in Theorem 
12.21 If /C < +oo, it is always satisfied by R = JC. If /C = +oo, there exists a 
sequence of invariant measures satisfying (}2"8")) with = max(A a , A&). Indeed, 
by Corollary EH M t N : = t" 1 f* /j N 'e sLn N^ds < u° R . Thus (M t N , t > 0) is 
tight, and one can pick any subsequential limit z/^ as t — > oo. In particular, 
( 128]) is satisfied if the invariant measure is unique for A r large enough, which 
holds regardless of /C in dimension one. 

2.2 Examples 

In the following examples, fl^ — {1> • • • , N — 1}, £7 = (0, 1), dQ = {0, 1}, 
Ajv(-) = A a lzn(-oo,i) zn(JV— l, +oo) > A(.) = A a l {0 } + A 6 l { i } , where A a , X b G 72. 
are given reservoir densities. For r\ G E, x, y G Z, a; 7^ y: o^r? is ?y — + S y if 
77(2;) > and rj(y) < JC, 77 otherwise; <5+?7 is 77 + ^ is 77(0;) < JC, 77 otherwise; 
#~77 is rj — S x if r)(x) > 0, 77 otherwise; cx^r/ exchanges contents of sites 2 
and y. If 7 G {5^, 0"^, cr^}, operator T G {A^, A±, E^} is defined by 



Misanthrope's process ( |15| ). Let p(.) be a probability measure on Z d , 
with finite support, satisfying the irreducibility assumption J2 n> i[p* n (-) + 
P* n (~ ■)] > 0; &(•) •) De a bounded function 01a E x E such that 6(0, .) = 0, 
and either (if JC = +00) b(n,rn) > for n > 0, or (if /C < +00) b(.,JC) = 
and 6(n, m) > for n > 0, m < /C. Then 



r/ = /o 7 -/. 




(29) 
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Attractiveness holds if 6(n, m) is nondecreasing in n and nonincreasing in m. 
Algebraic conditions 

b(n,m) b(n, 0)6(1, m) 

b(m + l,n- 1) 6(m + 1,0)6(1, n- 1) ^ 

b(n,m) — b(m,n) = b(n, 0) — 6(m, 0) 
ensure product invariant measures za, with explicit marginals. (fT5|) yields 



ff(p), #(p) := is p [b( V (x),r](ym Vx^V (31) 



m m 
N-2 

T J_T? r (\ ^ A+ _L TTt^°( AT _ 1 ^ \.U 

v jV-l 



The model b(n,m) = /(n)(l — af(m)) contains simple exclusion (/C = 1) 
and zero-range processes (/C = +oo, a = 0) , where / is any nondecreasing 
function on E such that /(0) = < /(l), a = l//(/C) if /C < +oo, or 
< a < l//(+oo) if /C = +oo. See jl] for other examples. Let 

6 + (p, n) = 9 p (m)b(m, n), 6 (n, p) = 9 p {m)b(n, m) (32) 
Then for p(l) = 1, 

x=l 

where the first term is the bulk generator (i.e. the part that does not in- 
volve interaction with reservoirs) and the next terms are the boundary part. 
For simple exclusion fl32l yields the usual entrance/exit rates ([351 H6| [16]: 
b + (p,n) = p(l — n), 6 (n, p) = n(l — p). For the zero-range process H in 
f[3"Tj) is increasing, thus the BLN boundary condition (|T51) - (j2Til) reduces to a 
Dirichlet condition at x — 0, and is void at x — 1. Microscopic counterparts 
are (since 6(n, m) depends only on n) the entrance rate independent of the 
particle configuration and the exit rate independent of the reservoir density. 

Asymmetric exclusion with overtaking (|3J). Let /C = 1, (e±, . . . ,e<i) 
be the canonical basis of R d , V = {±e 1 , . . . , ±e d }, k e N, (3J > for a e V 
and j = 1,. . . ,k, {I? + /3f a > (irreducibility), f3? +1 < (3? (attractiveness). 
Then 

k 

L = J2J2 c Uv)^o ja (33) 

aev j=l 
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where 

i-i 

<?(v) '= PjI[v(™) (34) 

n=0 

In words, a particle jumps to the first vacant site in a random direction 
a = ±ej with rate ft* if it is j < k sites ahead. u p is the product Bernoulli 
measure with mean p £ [0,1]. Thus f[T5]) yields 

h(p) = ^2h i (p)e i with hiW^pil-plYtifr-ft-j)?- 1 (35) 
t=i j=i 

By a proper choice of parameters, (I35p can produce any polynomial vanishing 
at p = and p = 1. Assume c? = 1, 2? = {±1}, k = 2, ft® = for a = —1, 
/3j = Then 

jV-3 



Ltt N) \ N (.) — 2_j L x + ftl^N-2,N-l 



x=l 
v2 



+ (ftx\ a + ft 2 K)At+ft 2 \ a r)°(l)Al 

+ [A(l - A 6 ) + ft 2 X b (l - A 6 )]A-_ X + ft 2 (l - A b )r/°(iV - 1)A N _ 2 
where the first line is the bulk generator and the sequel is the boundary part. 

Compound exclusion process. This model in the spirit of [12] illustrates 
multiple simultaneous jumps. Here d = 1 and /C = 1, 

L f = f o a 2i 3 o 5 ' 1 -f = (£ 2 ,3 o A 0)1 + S 2j3 + Ao,i)/ (36) 

which means that simultaneously a jump — > 1 is attempted and an exchange 
2^3 performed. By standard coupling methods of [H] one shows flTTj) . 
where u p is the product Bernoulli measure with mean p. Hence f fl5|) yields 
h(p) = p(l - p). Let Q N = {1, . . . , JV - 1}. Then 

AT-4 

Ln N ,\ N (.) = 2_s + Si i2 + Ajv_ 2 ,jv-i 
+ A a A|_ + (1 - A a )A* 

+ [1 - A.(l - 77°(1))]S 2 ' 3 + A Q (1 - ry°(l))(A+ o £ 2 < 3 + A+ + E 2>3 ] 

+ (1 - A 6 )A^_ 1 + (1 - \ b )r]°(N - l)(A A r_ 3 ,iv-2A 7V _ 1 + A 7V - 3 ,iv-2 + A 7V _ 1 ) 



+ A 6 (l - r)°{N - l))(A A r_3 iW _ 2 A+_ 1 + A JV _ 3iiV _ 2 + A 



iV-l) 
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where the first line is the bulk generator, lines 2-3 the left boundary part, 
and the sequel is the right boundary part. 



3 Proof of Theorem 12.1 

Let us define the microscopic 

entropy : ^ (rj(x) - ^(x))^ (37) 

x& d : \x\<r 

entropy production : D ± := LqJi (38) 



entropy flux : ^(77, £) := L 



x{t](x) - ^(x)) ± 

rrGZ d : |x|<r 



(39) 



where Lo = L is a monotone coupling of Lq with itself defined in ( 165|) . In 
Lemma 13.11 we show that D^ 1 < 0. The notational convention in (|3T1) -( |3"9~1) 
and herafter is that ± is replaced by + everywhere or — everywhere. In 
Subsection 13.11 we define a coupling Lq n ^ 1 (.),a 2 (.) °f two open systems in 
Qn with microscopic reservoir profiles A]y(.) and A^(.). We write Lq n X i n c 
when X%(.) = c. With this coupling we derive a microscopic counterpart to 
(TTTj). For < ip e C°°(M a! ) with compact support, let 

F±(V°,0 ■= N~ d ]T ^(x/N)( V °(x)-C(x)) ± (40) 



We show in Subsection 13.21 that 

+ N~ d dMx/^.TjfW^ ) 

+ C.N 1 -' Y ¥^(x/iV)(A] v (x)-A^(x)) ± 

+ C 2 N- 1 

(41) 

(see notations ©-Q), where C\ depends only on L , and C2 on L , fi, and 
sup {A^v (a;) : N E N, x & fl^}. (T4T1) is a generalization of the microscopic 



12 



entropy inequalities of [55], where the boundary term is the key novelty. 

Step two: mesoscopic entropy inequalities. Let be a coupled pro- 

cess with generator Lq n> x n (\ c and initial distribution jl No ' c 



defined in (ESJ), that is a coupling of fi No and u°; this particular choice of 
initial distribution will be used only in (|53]) ). Let < <p G C°°([0, +00) x R d ) 
with support contained in [0, +00) x K for some compact K C For 
i G [0, +00), let yjj = <p(i, .), and consider the mean zero martingale 



d s + NL UN , XNQ AF^(r ] T s ,^ s)ds 



In Subsection I3.3I we prove that, for every T > 0, 



?± ( No /•A ^ o^ 

(42) 



E 



(M^iip)) 2 <C(T,Sl,p,Lo)N- 



(43) 



Define 

K N± (<p) 



+ 



+ 



+00 



N~ d 9Mt, x/N)( v »°(x) - ffiWfdt (44) 
N~ d dMt,x/m-Tx]Hvm(x),&))dt (45) 
CxN 1 ^ (f(t,x/N)(f>f{X N (x))dt (46) 



(47) 



and /Cf ± (<p) = dSD + (USD, /C^) = (USD, /C^(^) = <m- By flU), (03 
(SB and D ± < 0, 



lim W,K, N± {<p)- = 

JV — »oo 



(48) 



for every e > 0. For i e {1,2, 3}, let JCf(cp : p) denote the i-th integral in the 
l.h.s. /C ± (v9,p) of (ITT)) , where (ft = <f>f. We show below that 



where 



lim lim sup E \)Cf ± (tp) - /Cf(y>, p N ' 1 ) 



p^(t, *) = AWtiS; ■= (21 + I)" 1 £ 



(49) 
(50) 
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This implies a mesoscopic version of (fT7j) : 

lim limsup EJC(ip, p 

For i = 2, fj49|) follows from (TTXj) . since 



(51) 



lim sup flUD < M / ^(t,x)0±(A(x))dtdH d - i (x) (52) 

N->oo J(0,+oo)xdQ 



where M depends only on L . For i = 3, (pf9j) follows from (1681) . since 

p 7 ^- lim flED = / ^(O,x)0±(p o (x))^ 
Jn 



(53) 



where /i^ ^ — lim denotes convergence in /i 7Vo ' c -probability. For (|49j) with 
i = 1, we prove in subsection 13 .41 that, for every e > and compact C M d , 



lim lim sup E 



lim lim sup E 



BD - 



(0,+oo)x< /iV 



%>(t, x)^(A [7Vx1 '^ )^^ 



(BSD - / d x <p{t,x)-i>t(A [Nx] ' l Vm)dtdx 

(0,+oo)x^ /JV 



(54) 



(55) 



Step three: passing to the limit. We skip here some technical details similar 
to [55]. A Young measure n(dt,dx,dp) ([TTj) is a nonnegative measure on 
(0, +oo) x M. d x 7?. whose projection on (0, +oo) x R d is the Lebesgue measure 
dtdx. It is called a measure- valued (mv) entropy solution to ffl6|) with data 
Po £ A G L°°(<9f2), iff. it satisfies the extension of fTTTj) . where the 

l.h.s. /C ± (<p,p) is replaced by (with = <pf) 



/C ± (v?,tt) 



[<9 t <p(t, x)0(p) + d x <f(t, x).ip(p)] n(td, tx, dp) (56) 

<p(0, x)(f)(p {x))dx 



M I ip(t,x)(f)(\(x))dtdH d - L (x 

'(0,+oo)xdQ 



p(., .) is an entropy solution iff. n(dt,dx,dp) = dtdx5 p ^^ is a mv entropy 
solution. For iV6N and Z e N, define the random Young measure 



it ,l (dt, dx, dp) := dtdxd p N,ifj- x ^{dp) 



(57) 
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By corollary \'S.11 (tt n,1 )m i is tight for the topology defined by convergence 
on functions F such that \F(t,x,p)\ < C(l + p)l^(t, x) for some C > 
and compact K. Since JC^fa, p N ' 1 ) = /C ± (<p, 71"^), any subsequential weak 
limit of 71"^'' as iV — > oo, / — - ► oo is supported on mv entropy solutions. By 
uniqueness of the latter ([61]). we have n N ' l (dt,dx,dp) — >dtdx5 p u x ){dp) in 
distribution, where p(., .) is the unique entropy solution. This implies, for 
every compact K C (0, +oo) x M. d , 

limlimsupE / \A [Nx] ' l r)%° - p{t,x)\ dtdx = (58) 

To conclude, we need to prove a^ t (dx) := a N (r]^)(dx) — > a t (dx) := pit, x)dx 
in probability for each t > 0. This follows from (l5T^j) if we prove 




lim lim sup P 



> S = (59) 



for every <p G C°°(Q) with compact support and every 5 > 0. (j5"9"|) is shown 
by considering martingale (H2j) c = and <pt = <p, so that 

The result then follows from (|4TI) . (H3|) and Doob's inequality. 
3.1 Coupling formalism 

We consider models for which L can be represented as follows: 

Lof(v) = [ lf(T u v) ~ f(v)]dm(u) (60) 



ju 

where (U, m) is a measurable space endowed with a finite nonnegative mea- 
sure m and, for each u G U, T u : E — > E has the following properties: it is 
local, i.e. there exists r G N (independent of u) such that T u r](x) = r](x) for 
all x G Z d \_B r (0), and T u r\ depends only on the restriction of r\ to Z D -B r (0); 
T" is conservative, moves a bounded number of particles and is monotone: 

E T >( x ) = E ^) ( 61 ) 

xeZ d : \x\<r xeZ d :\x\<r 

M := sup E \T u ri(x) - r](x)\ < +oo (62) 

r/GE, u£U „, . 
' z£Z d : |x|<r 

ry<^ T u r]<T u £ (63) 
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We set T x ' u = t x T u t_ x . Let us illustrate (16"UI) for models of Subsection 12.21 
Misanthrope's process. U = (0,sup6) x Z d , p(dvdz) = 1 

(o,supfe)^ ® p(dz). 

For (v,z) G U, T^ v > z 'rj = 5o !Z r] if v < b(r}(0),r)(z)); T^ v > z 'r] = r\ otherwise. 
Monotonicity assumptions on &(.,.) imply fl63|) . 

Asymmetric exclusion process with overtaking. U = (0, 1) x T>, p(dv, da) = 
dz®J2xev 8x{dcx). For a G V and 77 G E, let z a {rf) denote the first vacant site 
in direction a starting from (but excluding) the origin; set z a {rj) = if there 
is no such vacant site. For [v , a) G U, set: T^"' a 'r] = 5o :Za (ri)V if v < ft^r^', 
rp{v,a)^ _ ^ otherwise. (3f +x < ft® implies ([SSD- 

Compound exclusion process. U = {u}, p = 5 U and T u = cr 2 ' 3 o S 2 ,i- 

We couple n systems on Z d with common generator fT2l- fl60l) by 

£<») = ~ L (n) := ^(n^ (g4) 

L^firn,...,^) = [[f(T u r ]l ,...,T u r ]n )-f(r ]l ,...,r ]n )]dm(u) (65) 

Let x G K 1— > F p (x) := z/ p (?7(0) < x), with generalized inverse F~ (note that 
F~ l is nondecreasing wrt p by (1731) ). Let (U(x),x G Z d ) be i.i.d. uniform 
random variables on [0, 1]. Given S C Z d , define the probability measures 

<X),-,K) (F^(U(x)),...,F^ x) (U(x))) xeS (66) 



^(.),..^(.) ^m*)),---,^^*)))^ (67) 

where p 1 (.), . . . ,p k (-), resp. A x (.), . . . , A fc (.), are bounded 7?.-valued fields on 
S, resp. Z d \S. In ([6~4"|) . dHSD (EZD and below, we drop (n) when n = 2 and 
write c for A fc (.) = c or p fc (.) = c. Note that 

^(.),P 2 (.) " = (P 1 ^) - P 2 (^)) ± (68) 

for x G S, and similarly for the outer measure V. We define "mixed" coupled 
processes where some components may evolve in S, and others (indexed by 
0) in 7L d . For instance, 

^s!xH.),\H.)fif(v i,V°2,m) ■= J ^/W©^i^2©%^3)^5 3 ii ( . )iA 2 ( . } (^i,^ 2 ) (69) 
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is a coupling of L^iq, Ls,\ 2 (.) an d L. A subset of components is then 
Markovian with the generator defined from the sub-list of indexes. The 
process defined by f l69|) can be constructed from a Poisson measure uj on 
(0, +00) x S xUx (E zd \ s ) 2 with intensity 

£(dt, dx, du, dfji, df)2) = l[o t+oc )(t)dt © ls(x)dx®m(du) © i 7 ^^ \ A2 ( )(^7i, ^2) 

where (ix denotes the counting measure on Z d . For every (t, x, u, fj^ rj 2 ) £ k>, 
the configuration (77°, 7]%, 773) at time £— becomes at time i: 

fa?, %°, %) - ((T x ' u r]l © ry^ , (T x ' u r]2 © 7y 2 )° , T*%) (70) 
By the above construction, we have: 

Proposition 3.1 The subset of ordered configurations rji < ■ ■ ■ < rj n (or 
T)i < ■ • • < ffn) i s stable for the coupling defined by [64\)-{65}). or by ( TffPj) with 
A (•) < • - ■ < -^ n (-) containing no 0. 

Corollary 3.1 Le£ //° < (R > 0) be a probability measure on E s . Then 
Vt > 0, p°e tLs -^ < v° r , 

Attractiveness and irreducibility assumptions (see [UJ (U [151 E2] imply the 
following properties. Let I denote the set of invariant measures for L, and 
T the set of translation invariant probability measures on E 2 . Then there 
exist probability measures v p p i on E 2 for p, p' e 71, with marginals v p and 
Up>, such that (e denoting extremal elements) 

(jnf) e = (p, P ')e^ 2 } (71) 

p|p =► *P,A(v,Z)-- vttt = 1 ( 72 ) 

P < p' v p < v p , (73) 
3.2 Proof of ( gl] ) 

We first state properties of h ± , D ± , j and ? ± in (jT5]h (|57}, (|gg|l and fl39|) . 
Lemma 3.1 For every u EU and (r/, £) G E 2 ; 

/T(r/,0 = £+(£,77), ^(T^,T"0 < fr^O (74) 
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D ± < 0, D ± = on {i)<^}U{f)> f } (75) 
r(V,t)=3 + {t,V), J + fo,0 = l!L a(c\ lltl (76) 



i(^) - i(0 */ ^ > £ 

1^(17,01 ^Cmin^foO] (77) 

I^O-W.OI <C[%, 0] (78) 
/or a constant C = C(Lq), where h := h + + h~ . 

Proof of lemma \3.1[ (1751) — (1751) follow with elementary algebra from defi- 
nitions, (!60|) - (l63|) and (!7ij) . We prove the second statement in (IT3|) for 
which implies the same for /i~ by the obvious first statement. Let X- '■= V^^, 
and x+ '■= V V £. Since is nondecreasing, 

T u x- < T U V < T U X+ , T U X - < T u i < T U X+ 

Hence, 

£ (T^(x) - T«£(x)) + < £ (T»i7(x)-T" X -(x)) 

xeZ d :|x|<r xgZ d :|x|<r 

where the first equality follows from (I6TI) . □ 

Proo/ o/ gip. By dSSD , 

NL x F±{r), = NL x F± x {r), (79) 

where 

i&fo°,f ) == ^ E vfo/WG/) - f (v))* (80) 

For x G 

= N- d <p{xlN)T x h±{rf,?) +G± X ( V °,C) (81) 

where 

^,(^,0 = AT* £ {<p(v/N)-<p(x/N)){Ti (y) (82) 
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with 



2/e-B r (2) 

Hence (the expression below depends only on (77 , £°) because x G f2^) 



iVi^^f) = N 1 - d <p{z/N)T a D ± {r, ,?) 



+ NL x R±( V °,e) 



NL x F±( v ,0 = N 1 -*<p(x/N)T a D ± {ri,0 



N x - d ip{x/N)L a 



E fo(y)-e(j/)) = 

y£-B r (a;)\r2jv 



(83) 
(84) 
(85) 



produce the first two lines on the r.h.s. of (14"TI) . By (155]) and (16"2j) . 
contributes to the last line of (14"TI) . For x G <9f2^, by (1T9"j) - 



(86) 
(87) 



iV 1 -^^ £ [<p{v/N)-<p{x/N)Kviy)-Z(y)Y 



By Lemma EH (fSSD > 0. By (JMJ) and 



-L 9 



E Mj/)-£(j/)) = 



yeB r (x)\n N 



(88) 



<2Nm{U) E Mf) -£(*/)) 



± 



yeB r (x)\n N 



By integration w.r.t. z-Y^a 1 (.),a 2 (.) an d (JHSD , the contribution of fl8T|) is 
bounded by line 3 of (HIT) , plus an error that contributes to the last line 
of flUfl. By (J65D, ([HHD again contributes to the last line of (JH]). 



3.3 Proof of ( 1431 ) 

We first derive a bound on entropy production. 
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Lemma 3.2 For every < T < +00, 



sup E J f N l - d <p(t,x/N)T a (-D±)(Ti$i,t$)dt 1 < +00 (89) 



Proof of lemma Since M t Ar± (<^) has mean zero 



EF±(^,e^)-EF±(^ ,0=E \ T d t F^lU N m)dt (90) 

^0 



-IE 



Ln N ,\ N (.), c Fp t (rj N t, £ N %) 



tit 



(91) 



We bound ((9DJ) by Corollary O and apply (gTJ and Lemma O to (JSTJ. □ 
The martingale defined by (j42p has quadratic variation 



<M^)> 



N ^n w ,A w (.),c(^) - 2F^L QjViAjv( . )iC i v 



(92) 

Then, ( 1431) follows from Lemma [3.21 and Lemma [3.31 below. 



Lemma 3.3 For every nonnegative ip G C°°(R ) compact support, 
N |L njViAjv( . )jC (F±) 2 - 2F±L njViAjv( . )jC F± j (77 , £°) 



(93) 



Proof of lemmaWIE By (lM])-(j55jh 
Z(F±) 2 -2F±LF± = £ / [F^(T^,T^0-^(^0] 2 ^(«) (94) 



(in fact iGfi^U dfi^) with F± x defined in (ISO]). For x G dft^, by ([62]), 

[F± x (T^,T^0 - i^fa.fl] 2 < 4AT 2 MlN- 2d l K (x/N) (95) 
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for some compact K C M. d . These terms in (1941) contribute to iV d in (I93p . 
For x G by ([HI]), (JS2D and ([71]), the l.h.s. of (J2SD is bounded by 

-ANN~ 2d ^(x/N) 2 [^^(T^'V^'T) -^(v^C) 
+ W 2 N- 2d ~ 2 \\d x ip\\ 2 ^K(x/N) 

whose first term produces the first term in (1931) . and second term again 
contributes to N~ d in ([93]). □ 



3.4 Proof of (]54ft-(l55l) 

For every measurable subset K of [0, +oo) x Q, we define 



i^> N > i (v No ,e°) 



jK ,N \v No ,e o " 



K 



yeBi([Nx]) 



dtdx 



K 



(2i + iy d Yl ^ ± (^*(y),ejS(y))-^ ± (^'M?) 

We prove that, for every compact K and e > 0, with 7^ being X or ,7, 



dtdx 



lim hm sup EH™ (7T ,r°) = 



(96) 



By finite coverings we are reduced to K — (t — 5, t + 5) x -65(0:0) with 

JVo £No 



B s (x ) C f2. We consider the coupled process (r]N(t -s)+v ^N(t -s)+v '>£f') 



; (4) 



with generator L njViAjvQjCi0)0 



, see (ESD, and 77^ = © 0^, ^° 



~ £0^ extended to Z d so that ^ ~ v c . We claim that 



lim limsupEHf 5; xBi(xo) (r/:-,e: JV ) = 



I AT J 



(97) 



P|) follows from pTjl. (1751) and Lemma £00 below for (77^°, 77^') and (^ JVo , f 
We sketch the proof of the claim, similar to [55]. Write 



Ea/jv <| (2/ + l)- d £ - (A°' l V )) 

yeB t (o) 

(98) 
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M N = \Kl x [ T [Nx] p N e Nti dtdx (99) 

where p d ^ (r]' N ,£ N ), g±{ V ,t) = (17(0) - £(0))+ and Gf = <f>t in 

= j ± and Gf = ip^ 1 in J' ±,N, K has ^-marginal z/ c . By Corollary 

13.11 (Mjv)jv is tight for the topology defined by convergence on functions F 
such that \F(r), £)| < C(l + X^ert 7 ?! 2 -) + f° r some C > and finite 

T C Z d . By (1991) . every subsequential weak limit M lies in X fl T, hence 
by (J7T1) is a mixture of measures £ p y with imposed p' = c. limsupju . =0 (TUST) 
is the M-expectation of the quantity in braces. It is easy to see that the 
z> PjC -expectation of this quantity vanishes as I — ► 00 by (J75|) . (1751) . (1T5|) . (Tl8l) 
and the spatial ergodic theorem. 

Lemma 3.4 Let B$(xo) C f2, and (x^°, C/O a coupled process with generator 
Lfijy,A w (.),0 anc? < z/£, ; < z/r m distribution for some R E 71. Then, 
for every t > and < 5' < 5 - Vt, with V = V(L ) and C = C{6, 5', L ): 

VN~ d Yl \x%°(x)-a(x)\<TEN- d J2 \x% o (x)-( N (x)\+CN^ 

x€NB s ,(x ) x£NB s (x ) 

Proof of lemma \3Jh Let V > 0, e > 0, and (with the notation of ([40]) ) 



Me := \\ x \ + £ > (fi(t,x) := H * 5 £ (5 — Vt — \x — x \ e ) 

F(t, X ,C):=F+ ti . ) ( X ,C) + F- M ( X ,C) 

where H is the Heaviside function and 5 £ a standard mollifier. Using ()42l) for 
F, dHJ) and (E7J), we have 

EF(t,£$,C&) < EF(U ^ C) + CAT- 1 

+ E r^- d ^{^,x) + C'|^( S ,a;)||^ s -C(x)|}rf S 

for some constant C > 0, where expectations exist by Corollary 13.11 Since 

dttp(t,x) + V\d x tp(t,x)\ < 
we obtain for large V (which implies the result for small e) 
EF(t, Xm, Cm) < EF(t, , Co") + CiV" 1 * 

□ 
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4 Proof of Theorem 12.3 



Proof of (i). ( 128]) implies tightness of (a N (rf)') „. Given p E 1Z, we set 



zy° : = ^e^n^.A^t.). By ((2SD, Proposition O and stationarity of v No , 



<°m < < (100) 

Let M No := v No o Mj° := z/J o a" 1 . We consider on the set of Radon 
measures on Q the partial order defined by a < (3 iff. /3 — a is a nonnegative 
measure. Since r]° i— > a N (i]°) is nondecreasing, 

M % < M^ < Mr° n1 . (101) 

for the associated stochastic order. By Theorem 12.11 any subsequential weak 
limit M° of M^ as iV -> oo satisfies 

S p o(t,.)dx <M°< 5 p R( tt ,) dx (102) 

for every £ > 0, where p r (., .) denotes the entropy solution to (flBT) with data 
Po(-) = r in O and A(.) on <9f2. By Theorem I2.2[ p°(t, .) and p R (t, .) converge 
to iift(.). n (A a , A&) as t — > oo in L 1 1 oc (f2" 6 ). £ — > oo in (j!02p shows that the 
restriction of M° to fi£ 6 is 5 fih( n ( Aa ,A b )^- 

Proo/ o/ (it/ Set c := %.). n (A a , A 6 ). Let (r/f °, v $, rffi, £f °) be coupled 
with generator £^ iAjv( . )iAjv{ . )iAjv( . )iC , V°,o ~ < f o r r e {0,^}, ^° ~ and 

0^ = <o < rff° < nSfi ( 103 ) 

a.s., which is possible by (175]) and (Eg) . We apply (JH]) to (r^°, and §HD 
to r]^° (here with limiting Young measure dtdx5 p r( t ^(dp)). This yields 

jimE(/ T / (f 7p %([iVx])-^([^])) ± &:*)= T / (p r (t,^)-c) + ^t 



JK J JO JK 



> 



for any T > and compact K C f2™ 6 . By Theorem 12.21 applied to p r 

lim lim E {t- 1 / T / (v?M[Nx)) - dxd£ = o) (104) 

T->ooN-+oc I Jo JK J 
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for r G {0,R}. By ( TTU3D and Proposition EH < Vm < Vrnv Thus 



lim lim E < T -1 



T 



-/if 



^ ([iVx])-^ ([iVx]) \dxdt = 



(105) 



Since g is bounded and local, \g(rf) — g{£)\ < C J2 x eB,(o) \v( x ) ~ f° r 

imply, for some 



-No 



xeB t (o) 

y No 



some C > 0, I G N. Stationarity of z/^ and £ 
compact if' C (which concludes the proof by (11051) ). 



\g N (x) - g\ dx 



T -i 



K 



K 



E 



< CE^T" 1 



[t\Nx]9 {Vm) ~ T [Nx]9 (Zm)] dt 
\Vm([Nx])-^Nt([Nx])\dxdt 



JK 



5 Proof of Theorem 2.2 



We need comparison results for entropy sub/super solutions. An entropy sub- 
solution, resp. super-solution ( [10J ) to (fTEl) . is defined by restricting (fT7|) to 
{<j>,1>) = (0+ resp. (0,^) = (0-,O, and lies in C°([0, +00); L} oc (Q,)) 
by [52]. 

Proposition 5.1 Lei p 1 (., .) (resp. p 2 (., .)) be an entropy sub-solution (resp. 
super- solution) to ( 175]) in /or date Po(.), Po(-) ^ n ^ anc ^ ^i(-) — ^2(0 i n 
<9f2. Tnen ; setting V = max (| |/t' o pj| 1^ , | \h! o Pol loo)' 



(p 1 ^) - p 2 (t,x)) + dx < / (pj(ar) - p 2 (x)) + dx (106) 

sinB(i ,R) JnnB(x ,R+vt) 

for every t > 0, Xq G Q and R > 0. In particular (maximum principle), 
/(.,.) < max(pj(.)) VmaxAiO), p 2 (-, •) > min(pg(.)) A min A 2 (.) (107) 

s2 OS! i2 ail 

(11061) follows from Lemma 2 of [M] by standard arguments, and (1 1 7j) by 
taking constants as particular entropy solutions. The next one is proved in 
Appendix IA.U 

Lemma 5.1 Let p(., .) be the entropy solution to [To]) in Q± C K d itfit/i data 
Po(.) in Oi, A(.) on c?fii. Lei ^2 ^ e an open subset of £ C <9fii fl (9^2? 
and i? constants such that R~ < p < R + a.e. in 

(0,+oo)xfi 2 - Let A ± (.) 5e 
egna/ to A(.) on £ and R^ on <9fi 2 \X. Tnen p(., .) is an entropy sub-solution 
(resp. super-solution) to [W) in f2 2 for data po(-) restricted to Q2 and A + (.) 
(resp. } ~ (■)) on 00,2 ■ 
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5.1 The case Q = Vt n a b 

By Proposition 15. 1[ the entropy solution p(., .) lies between entropy solutions 
with the same boundary datum and uniform initial data and | |po| loo- Thus 
it is enough to consider uniform data p (.) = r G 1Z on fi™ 6 . Let p n (t,y) be 
the entropy solution to 

d t p{t,y)+d y h n {p{t,y)) = (108) 

on (a, 6), with /i n (.) := h(.).n, data po(.) = t in (a, 6), and A (7) = A 7 on 
{a, b}. It is easy to see that 

Lemma 5.2 p(t,x) := p n (t,n.x), defined for (t,x) G (0, +00) x VL™ h , is the 
entropy solution to [To]) with data po(.) = r in and [2^) on dQ™ b . 

We are now reduced to proving 

Proposition 5.2 Let f G C 1 (7?.) ; and p(t,x) be the entropy solution to 

d t p + d x f(p) = (109) 

on (a,b), with data po(.) = r E 1Z in (a,b) and A(7) := A 7 on {a, b}. If 
Mf(\ a , A&) = {Rf(\ a , At)}, then p(t, .) — > R f (X a , A&) m L}{{a, b)) as t -> 00. 

The following particular case will be proved in Subsection 15.21 

Lemma 5.3 In i/ie context of Proposition \5.Sl assume X a = At = A with 
/'(A) 7^ 0. T/ien, for any r ElZ, p(t, .) — > A in L 1 ((a, 6)) as t — > 00. 

Proof of proposition 1 5. B We treat A a < A&, the other case is similar. For 7 G 
{a, 6}, t > 0, x G R, let p 7 (t, x) be the minimizer of /(p) — p[x — (a + 6 — j)]/t 
on J 7 , with J a = [A{, A&] and = [A a , A{]. We claim that the restriction of p 7 
to (0, +00) x (a, b) is the entropy solution to (11091) with uniform initial datum 
A^ in (a, 6), and boundary datum A(.) on {a, 6}. Indeed: one the one hand, 
by [3 Proposition 4.2], p 7 is defined (i.e. the minimizer unique) a.e. and is 
the entropy solution to (11091) on R with initial datum min(J 7 )l(_ 0Oia+6 _ 7 ) + 
max(/ 7 )l( a+b _ 7 +0O ), whose restriction to (a, 6) is A^; on the other hand, the 

f( 1 

definition of p 7 (., .) easily implies BLN boundary conditions p 7 (t,a + ) G £ x 
and p 7 (t,o~) G £\ ' with the sets defined in (]2"U]). This proves the claim. 
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Since Rf(X a ,Xb) is the unique minimizer of / on either interval 7 7 for 7 G 
{a, b}, we have 

lim p 7 (t,x) = R f (X a ,X b ), \/xe(a,b) (110) 

By definition of A{ and X{, there are sequences (A*), (A£) such that 

A* T A{, A fe fe | A^, f'(X k a )^0, f'(X k b )^0 (111) 
as k —>■ 00. By Proposition 15.11 and semigroup property, for < s < t, 

[ (p{t,x) - p\s,x)Y~<dx < [ (p(t- s,x) - X^dx (112) 

J{a,b) J(a,b) 

with signs a a — — , a b — +. For 7 G {a, b}, let pj.(t, x) be the entropy solution 
to fllQ9p with uniform data r in (a, b) and A^ on {a, b}. By Proposition 15.11 

p a k (t-s,.)<p(t-s,.)<p b k (t-s,.) (113) 
By Lemma [5T31 applied to p^, 

lim p^(£ - s,x) = A^ inL 1 ^)) (114) 

We take lim sup s ^ lim sup^ in (jTT2l> . By ffTTOl for the l.h.s, ffTT3|) -ffTT4l) 
for the r.h.s., 

limsup f {p{t,x)-R f {X a ,X b ))^dx<{b-a){X^-X! / )^ 

t->oo J(a,b) 

and we conclude with k — > 00 and (II lip . □ 
5.2 Proof of Lemma I5J3] 

The result is trivial for r = A: then p(t, .) = A. We treat r > X, r < X being 
similar. Let p(t, x) be the entropy solution to (11091) on IR with datum 

p (x) = Al(_ 00ia ) U ( 6i+00 )(x) + rl (a>6 )(x) (115) 

p (.) has bounded space variation, thus also p(t, .) for t > (see e.g. [57] or 
|63j). Hence, limits p^x 1 * 1 ) exist. By maximum principle (I107p . 

p(t,.)>A, p(t,.)>A (116) 
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for every t > 0. By (11161) and Lemma 15.11 applied to p, with Qi = R, 
fi 2 = (cl, b), S = and _R~ = A, the restriction of p to (a, 6) is an entropy 
super-solution to (I109p in (a, b) for the same initial and boundary data as 
p. Thus, by Proposition I5.1[ p(, .) < p(t, .) for every t > 0. To conclude we 
show that, for any bounded interval JcK, 

p < A + e in I for large enough t (H7) 

We prove (11171) as follows. The initial condition (11151) crosses X + e at a and b. 
In Proposition IA.ll we prove that there are uniformly Lipschitz trajectories 
< £2(0 such that xi(0) = a, £2(0) = b, and for all t > 0: p(t, .) < A + e 
on (—00, xi(t)) U (x2(t), +00), p{t, .) > A + e on (xi(t), X2(t)), and xi(t) = 
£2 (t) implies Xi(s) = X2{s) for all s > t. The latter implies (11171) . We claim 
that Xi(t) < x 2 (t) for all t > implies lim i ^ 00 xi(t) = +00 if f'{\) > 0, or 
lim^oo X2{t) = —00 if /'(A) < 0, whereby All Tj) follows. We prove the first 
part of the claim, the second one being similar. Set 

_.„ (A , A+£): . mf nx + *)-m (118) 

v ' ee[x,x+e) X + e -9 v ' 

v(X, X + e) > for small e. Let 5 > such that 

5e> I (p (x) - (A + = (b - o)(r - (A + e)) (119) 

Set x(t) := a — 5 + vt, and 

T := sup{r > : x{t) < Xt(t) for every t G (0, r)} (120) 
We prove that T = +00. Assume T < +00, and let R > b such that 

x 2 (t)<R, Vte(0,T) (121) 
which is possible by Lipschitz continuity of a^O- Set 

I(t)= (p(t,x) - (X + £))~dx, J(t)= {p(t,x)-(X + e))~dx 

Jx(t) Jx 2 {i) 

The measures defined on (0, +00) x R by 

m ± (dt, dx) := d t <j)± +e (p(t, x)) + d x ip± +£ (p(t, x)) (122) 
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where ((f>f,i/jf) is given by (TTHT) with h = f, satisfy m + = m < 0, by (11091) 
and (HID with Q = R. By integration of (jT211 . 

7(T)-I(0) = m{(t,x) G (0,T) xl: xE [x(t),R]} 

+ J [ip^ +£ (p{t,x{t)-))-v^ +£ (p(t,x(t)-))]dt ^ 123 ^ 

^ +£ (p(t,i? + ))rft 







(124) 



J(T)-J(0) = m{(t,x) G (0,T) x R : a: G [x 2 (t),R]} 

+ [ fc e (p(t,x 2 (t)"))-x 2 (t)0^ £ (p(t,x 2 (t)-))]^ 

Jo 

By (H20D, 

p(t,x(t) ± ) < A + £ (125) 
for t G [0,T). dUSD, (HU) and ([125]) imply 

^(fi{t,x(t)-))-v<f>^(ji(t,x(t)-)) > (126) 
By definition of ^(t), 

p(t,x 2 {t)~) > A + £ 

which implies 

^ A - +£ (p(t,x 2 (t)-))-x 2 (t)0^ £ (p(t,a; 2 (t)-)) =0 (127) 
Since x(T) = x\(T) by (11201) and continuity of x(.) and Xx(.), 

pR pa pR 

I(T)-I(0) = / (p(t,x)- \-e)~dx- ipo(x) - A — e)~dx — I (p (x) - \ - e)~ dx 



Jx 2 (T) Ja-8 Jb 

= J(T) - J(0) - / (p {x) - (A + e))~dx < —A — 5e 

J a— 8 

(128) 

where A := f Q i[)^ +£ (p(t, R + ))dt, and the inequality in ( 1128ft follows from 
(1121) . (fT27l) and m < 0. Let 

p rb+Vt 

K(t) := (p(t,x) - (X + e)) + dx = / (p(t, x) - (A + e)) + rfx 



a-Vt 
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where the second equality follows from Proposition 15.11 which implies 

p{t,x) = X, Vx eR\{a-Vt,b + Vt) (129) 

Integration of (I122p combined with (11291) yields 

K[T) — K(0) = m{(0, T) x E} (130) 

Using (Tim (jngj) . (IT231 . (pSj) and m < 0, we can contradict (Tim 

-A-fe < -A + K(T) -#(0) = -A + m{(0,T) x M} 

< -A + m {(t, x) e(0,T)xl: 16 [x(t), R)} < I(T) - 1(0) 

(131) 



5.3 The general case (1221) 



The following is an analogue of Lemma 15 .31 

Lemma 5.4 Let p(., .) be the entropy solution to [To]) with uniform data 
r G TZ in Q and A G 1Z on dQ, where h'(X).n ^ 0. Then p(t, .) — > A as 
t — > oo m L[ oc (0). 



Proof of lemma \5Jj\ Denote by p a ',b'(t,x) the entropy solution to (flE]) in 



^a' 6' w ith uniform data r in f2™, v and A on <9f2", fc ,. By Lemmas 15.21 and I5"T5} 

p a /, y (*,.)-A, in LLW^O (132) 

We consider the case r > A, the case r < A being similar. By maximum 
principle (I107p . 

P(t,.)>\, p a >,b>(t,.)>\ (133) 

for every t > 0. By (1 1 3 3 [) and Lemma I5TT1 with p = p a ',b', &i = ^"'fe'; 
^2 = ^, S = 0, R~ = A, the restriction of p a > ,b> to (0, +oo) x Q is an entropy 
super-solution to (I16p in Q for the uniform data r in Q, X on <9fi. Thus, by 
Proposition 15. 1[ p a ',b'(t, .) > p(t, .) in f2 for every t > 0. This, (11331) and (11321) 
imply the result. □ 

We now conclude the proof of Theorem 12.21 in the general case fl22l) . We 
consider A a < A&, the reverse case being similar. We set /(.) = h(.).n. We 
denote by T t the solution semigroup for ffTBT) in Q boundary datum A(.); thus 
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p(t,x) = T t po(x). Proposition 15.11 implies that T t is monotone. We consider 
sequences satisfying f illip and define, for 7 e {a, b} and k G N, the following 
entropy solutions to (JIB]) , p 7 has uniform initial datum A^ in Q, and bound- 
ary datum A(.) on dfl. p a has uniform initial datum A{ in and boundary 
datum A a (.) on dfl™ b , equal to A{ on {n.x = a}, X b on {n.x = b'}. p b has 
uniform initial datum X( in fi", 6 ; and boundary datum X b (.) on <9fi" b/ equal 

to A a on {n.x = a'}, X( on {n.x = b}. pi has initial datum p (.) in Q, and 
uniform boundary datum A^ on dfl. (11131) still holds in the present setting 
by Proposition 15.11 By Lemma [5.21 and Proposition 15.21 for 7 G {a, b}, 

lim p 7 (t, .) = R f (X a , X b ), in LUK, b ) (134) 

t— >+oo 

By Lemma 15.41 and Proposition 15. 1[ 

lim p2(t, .) = A? inLL^) (135) 

We claim that: (i) the restrictions of p a and p a to fin£7™ b / are respectively a 
super-solution and a sub-solution to (TTB"]) in n f2™ 6 / for the uniform initial 
datum A{, and boundary datum A a (.) equal to A{ on {n.x = a} and X b on 
dfl b ; (h) the restrictions of p b and p b to f2nf% b are respectively a sub-solution 
and a super-solution to (TTE1) in fl fi™, 6 for the uniform initial datum A{, 

and boundary datum A fe (.) equal to A a on dQ a and A{ on {n.x = b}. Indeed: 
first note that, by maximum principle (11071) . 

X{ < min(p a (t, .), p a (t, .)) < max(p a (t, .), p a (t, .)) < X b (136) 

for every t > on f2 fl Then apply Lemma [5. 11 (a) to p, with Qi = Q, 
Q 2 = n fi^, S = 9^6, 9fi 2 \S = {n.x = a}, R~ = X f a , (b) to p a , 6 / with 
fti = d 2 = n fi^, S = {n.x = a}, <9ft 2 \£ = <9ft 6 , i?+ = X b . This 

establishes (i),(ii) being similar. By Proposition 15.11 for 7 e {a, 6}, 

(p(i,x)-p 7 ( S ,x))^ < / (T s p*(t- S ,x)-p 7 ( S ,x)r^:r 

< I {T s p k 1 {t-s 1 x)-T s X)y-tdx+ ! {T s X s 1 -p\s 1 x)Y-dx 



< 



n™ b nB s (x ) Jn* b nB s (x ) 
(p*(t- s,x) - X}Y-dx 

W b nB s (x +Vs) 

(137) 
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where < s < t, a a = — , crj, = +, Bs(x ) C Q. In (I137p . the first inequality 
follows from (11131) and monotonicity of T s . The first term on the second line 
is bounded by the third line by Proposition 15.11 The second term on the 
second line, where T S A^ = p 7 (s, .), is zero by Proposition 15.11 and the above 
claims (i), (ii). The result follows by lim^oo lim sup^^ lim sup^^ in (11371) 
and ffTMl)- flT35l) 



A Proofs of some PDE results 

A.l Proof of Lemma [5JJ 

We may assume in (fT7|) that 

M>sup{|fc'(r)|: 0<r< ||p||J- (138) 

Let ((fi ± ,ip ± ) = (4'ti' l l ) t)- By (jTTl) . a nonpositive measure is defined by 

m ± (dt, dx) := dt^ip) + V,.^ ± (p) 

on (0, +oo) x Qi. Let O be an open subset of Q± with a boundary of locally 
finite perimeter. By [T21 Q3], there is an inner normal trace V>c>(t, x ) on 
(0, +oo) x dO for ^(pit, x)), that satisfies Green's formula 

Io(<p) ■= / {^{pit, x))dt<p(t, x) + ^(p(t, x)).V x (p(t, x)]dtdx 

ip{t,x)m ± (dt,dx) - I ip{t,x)^{t,x)dtdn d ~ l (x) 

0,+oo)xO J(0,+oo)xdO 

- / ip{0,x)<f) ± (p {x))dx 
Jo 

(139) 

for every smooth tp > with compact support in [0, +oo) x W d . For O = Q\, 
with (fT7|) and y? supported in (0, +oo) x M. d , we obtain 

^ i (t,x)<M0 ± (A(x)) (140) 
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Then, with O = Q 2 in iTB . 

In 2 (<P)+ / <p(0,x)(f) ± (p (x))dx + M ip^x^C^ix^dtdH 
Jo.2 i(o,+oo)xaa 2 



> 



+ 



(0,+oo)x£ 



M0 ± (A(x))-^ 2 (t,x 



n 2 < 



dtdH^Ux) 



(0,+oo)x(dQ 2 -E) 



M(j) (R ) — IpQ (t,x) 



dt dTi 



d-l, 



x) 



( 141 ) 

The second line of (11411) is nonnnegative because of (I140p . since $q = ^ 
a.e. on (0, +00) x (dOi fi <9f2 2 ). For a.e. (t,x) G (0, +00) x £7 2 and any 
unitary vector n G R , i?~ < p(t, x) < R + and (I138p imply 

^(p(t, a:)). 71 < Af0 ± (p(t, x)) < M0 ± (i? ± ) (142) 

Thus ^ 2 (t,x) < M0 ± (i? ± ) a.e. on (0, +00) x <9f2 2 , whence the result. 

A. 2 Sign changes for entropy solutions 



Proposition A.l Let c G 7Z, and p be an entropy solution to U09\) on 
with initial datum po(.). Assume there exist 



00 = x < Xi < ■ ■ ■ < x n < +00 = X n+ i 



(143) 



such that on (xj__i,£j), Po < c for odd i and p > c for even i. Then there 
exist trajectories t > 1— ► Xi(t) for % = 0, . . . , n + 1 such that: 



(i) -00 = x (t) < x\(t) < ■ ■ • < x n (t) < x n+ i(t) = +00 

(ii) Xi(0) = Xi for i = 1, . . . ,n 

(Hi) Xi(.) is Lip schitz- continuous for i = 1, . . . , n. 

(iv) On (xj_i, x^ we have p(t, .) < c for odd % and pit, .) > c for even i. 

(v) If Xiit) = Xi + i(t) for some t > 0, then Xi(s) = Xi+i(s) for all s > t. 

We use Godunov's scheme (see e.g. [20]) to prove the result for n = 1, 
then extend it to larger n by finite propagation property. Let At > and 
Ax = e > satisfy the CFL condition (see |30j): 

A := At/ Ax < 1/(2V), V = swp{\h'(p)\ : < p < ||p |L V c} (144) 

For j e NU(N+ 1/2), define pj(.) = p|(.) as follows. p (.) is the exact initial 
datum. For j G N, on each interval (kAx, (k + l)Ax), pj + i/ 2 (.) is equal to 
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the mean value pj^ of pj(.) on this interval. For j G N, Pj + i is the entropy 
solution to (I109p at time At with initial datum Pj+i/2- By (11441) and Proposi- 
tion 15. 1| pj+i is a superposition of the noninteracting Cauchy problems with 
data Pj,k-il{-oa,kAx) + Pi,fcl(fcAx,+oo) (Riemann problems). 

Step one: n — 1. We construct a trajectory j i— > £/i(j) = yi{j), where 
j G NU(N+l/2), such that: (ii') yf(0) = a*; (hi') ^(j + 1/2) - < Az; 

(iv') pj(.) < con (— oo,?/i(j)) and p.,(.) > c on (yi(j), +oo). We then de- 
fine x^(t) = |/f([t/At]). By (hi') and (I144p . xf(.) converges uniformly to a 
Lipschitz-continuous £i(.). pf t /At](-) converges (see [30]) to p(t, .) in /^(IR) 
as £ —>■ 0. This and (iv') imply (iv). 

We now construct yi(j). We set yi(0) = x\. Let j G N and = [yi(j)/Aa;]. 
Assume ?/i(j) satisfies (iv'). We set yi(j + l/2) = kAx if p^ > c, = (/c + l)Ax 
if p^fc < c, or of these values if pj^ = c. Then (iii')-(iv') hold. Assume 
yi(j + 1/2) = kAx for some fceZ, and 

Pj,fc < c for j < k, pj^ > c for j > k (145) 

Let I G Z. By construction of the Riemann problem (see e.g. [30]), on 
((/ — l)Ax + VAt, (/ + l)Ax — VAt), Pj+i(.) is a monotone function, with 
constant value Pjj-i on ((/ — l)Ax + Vt, lAx — Vt) and constant value p^i 
on {I Ax +Vt, (I + l)Ax - Vt). This and ffT45l) imply that: 

(a) For / < k, we have p j+ i(.) < c on ((/ - l)Ax + VAt, (/ + l)Ax - VAt) 

(b) For / > fc, we have p i+ i(.) > c on ((/ - I) Ax + VAt, (I + I) Ax - VAt) 

(c) For I = k, there exists y G (kAx — Vt, kAx + Vt) such that Pj+i(.) < c 
on ((k - l)Ax + Vt, y) and p i+ i(.) > c on (y, (k + l)Ax - Vt) 

We set yi (j + 1) := y in (c) above, (a)-(c) and (11441) imply (iii')-(iv'). 

Step two. We prove the result for n > 1 on [0, T) with T > 0. For 
i G {1, . . . , n}, set 

Po — Polfo-i^i+i) + cl]R\( Xi _ 1)a;i+1 ) 

Let p*(t, x) be the entropy solution to (I109P with initial datum p l . p l satisfies 
the assumptions of the proposition for n = 1 with a single initial interface at 
Xi (modulo the symmetry x i— > —x for even i, which changes / in (11091) into 
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— /). For each i = 1, . . . ,n, we obtain a Lipschitz trajectory t > i— > 
such that (a) Xj(0) = x,; (b) for odd i, p*(.) < con (— oo,Xj(t)) and p l (.) > c 
on (x{(t), +00); (c) for even z, p*(.) > c on (— oo,Xj(t)) and p*(.) < c on 
(xi(t),+oo) By Proposition 15. 1[ 

p\t, .) = p[t, .) on (x^x + Vt, x i+1 - W) (146) 

By (b)-(c) above, f!146j) and Lipschitz continuity of £»(.), there exists T > 
such that Xj(t) < x i+ i(t) and (iv) of the proposition holds for t < T. 

Step three: conclusion. Let T* > be the supremum of r > such that 
the proposition holds on [0,r). Then it holds on [0, T*]. Suppose T* < +00. 
Partition {0, . . . , n + 1} into maximal subsets Iq, . . . , I m +i, where 1 < m < n, 
such that (i) max 1^ < min for each k — 0, . . . , m; (ii) Xj(T*) has constant 
value (denoted by y^) on each I k . Note that J = {0} and I m+ i = {n + 1}. 
Thus 

-00 = y < yi < . . . < y m < y m+1 = +00 

satisfy the assumptions of the proposition for initial datum p(T*, .). By step 
two, for some T** > 0, the result holds on [0, T**) with trajectories yj(.) 
(for i = 0, . . . , m + 1), initial sequence (yi, i = 0, . . . , m + 1) and initial 
datum p(T*,.). Let a^t) = y k (t) for i e I k for T* < t < T* + T**. By 
semigroup property the result still holds on [0, T* + T**) for the sequence 
(xi, i — 0, . . . , n + 1) and datum po(-) : a contradiction. 
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